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Abstract 

o; 

O . We present a complete and optimal set of observables for the exclusive 4-body 

B meson decay — )• K*^{-^ K7r)i^i~ in the low dilepton mass region, that 
contains a maximal number of clean observables. This basis of observables is built 
in a systematic way. We show that all the previously defined observables and any 
observable that one can construct, can be expressed as a function of this basis. This 
set of observables contains all the information that can be extracted from the angular 
distribution in the cleanest possible way. We provide explicit expressions for the full 
and the uniangular distributions in terms of this basis. The conclusions presented 
here can be easily extended to the large-g^ region. We study the sensitivity of the 
observables to right-handed currents and scalars. Finally, we present for the first 
time all the symmetries of the full distribution including massive terms and scalar 
contributions. 



1 Introduction 



The rare decay — K*^{-^ K7i)i~^i~ provides unique opportunities in the search for 
New Physics in flavor physics due to the wealth and variety of angular observables avail- 
able experimentally. A total of 12 different angular coefficients characterize its angular 
distribution, each being a function of the invariant squared mass of the lepton pair, g^. 
Although a full angular analysis with a small binning requires a good deal of statistics, 
it constitutes a conceivable goal for LHCb, at least in its upgraded form. 

First data on the decay rate and several angular observables are already available from 
BABAR, Belle, CDF and LHCb. BABAR [T] has measured the decay rate, the forward- 
backward asymmetry Ape and the K* longitudinal polarization fraction Fl, all integrated 
separately in the low and high-g^ regions. Concerning g^-dependent measurements, Belle 
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[2] has provided a measurement of the total branching ratio and Aps, while CDF [3] has 
provided also measurements of Fl as well as the observables Aj,' and A^^ (see [HIS]), with 
a measurement of the dependence in the form of 3 bins in the low-g^ region (below the 
J/ijj resonance), a bin between the J/ijj and ijj' resonances, and two bins in the high-g^ 
region (above the The LHCb collaboration has also provided measurements of the 
branching ratio, Ap^ and based on ~ 300 pb~^ of data [6], while a larger data set 
of 1 fb~^ is on tape. In order to cope with limited statistics, near future plans focus on 
fully integrated observables, where the dependence is lost, at least within the low- and 
high-g^ regions. 

On the theoretical side, the interest is focused on the tayloring of observables with 
desired properties. These properties are: 1) a reduced hadronic uncertainty, and 2) an 
enhanced sensitivity to short distance contributions from New Physics (e.g. right handed 
currents, etc). Concerning hadronic uncertainties, the objective is to minimize the de- 
pendence on the soft form factors, which are difficult to compute and are the source of 
large theoretical uncertainties. This is achieved with the construction of ratios of angular 
observables where a complete LO cancellation of the form factors occurs. The search for 
observables with such desired properties has led to the formulation of a set of observables 
called aP g], Ap^'^^ 011] and A^^"''""^ |8] at \ow-q^, and an analogous set H^^^'^'^'^^ [9] 
at high-g^. These observables have been studied in detail and they indeed exhibit a low 
theoretical uncertainty and a clean sensitivity to characteristic New Physics features^] 

The source of experimental input is the differential decay distribution of the 4-body 
final state A'*°(— K'K)i~^i~ . It is described by four independent kinematic variables, 
which are traditionally chosen to be: the invariant squared mass g^ of the lepton pair; the 
angle 6^ between the directions of fiight of the kaon and the B meson in the rest frame 
of the K*^; the angle 6i between the directions of fiight of the i~ and the B meson in 
the dilepton rest frame; and the azimutal angle between the two planes defined by the 
lepton pair and the Ktt system^. In terms of these kinematic variables, the differential 
decay distribution can be written as 



dq"^ dcos 6k dcos 6, 



9 
3'2^ 



Jis sin^ 9k + Jic cos^ 9k 



■ {J2s sin^ 9k + J2c cos^ 9k) cos 29i 



^3 sin^ 9k sin^ 9i cos 2</) + J4 sin 29k sin 29i cos + J5 sin 29k sin 9i cos 
{Jqs sin^ 9k + J&c cos^ 9k) cos 9i + J7 sin 29k sin 9i sin (p + Js sin 29k sin 29i 

Jq sin^ 9k sin^ 9i sin 20 



sm I 



fl) 



The explicit dependence of the coefficients Ji(g^) in terms of transversity amplitudes (Ai) 
is given in Section |2l The point to emphasize here is that only observables that respect 



^For a representative set of references discussing the phenomenology of this decay mode see [TUl fTTl 

[ia[Ta[Ti[i5i[ii[ni[ii[i9]. 

^This definition of the kinematic variables coincides exactly with that of Refs. [HI HD] 



2 



the symmetries of this angular distribution can be obtained. These symmetrieqj are 
transformations among the transversity amphtudes that leave invariant the coefficients 
Jj(g^) of the angular distribution [7]. The number of amplitudes and the number of such 
symmetries determine precisely the number of degrees of freedom that are available from 
an angular analysis alone. This number is given by [7] 

nexp = 2 - , (2) 

where Uexp denotes the number of experimental degrees of freedom, ua is the number of 
transversity amplitudes Aj (j = 1, . . . ,71.^) and is the number of continuous transfor- 
mations (or symmetries) of the A's that leave the J's invariant. 

The number Ug of continuous symmetries can be obtained by inspection of the set 
of infinitesimal transformations of the A's, by the method described in Ref. [7]. In this 
way, one can infer the true number of independent experimental degrees of freedom nexp- 
If the number rij of coefficients Jj is larger than riexp, then the Jj are not independent 
observables, and a set of (nj — nexp) relations must exist between them. As shown in 
Ref. [7], there are exactly 8 independent J's (and consequently 8 independent observables) 
in the case of massless leptons, and 10 considering the mass terms -not including the 
CP-conjugated mode-. Adding scalar contributions increases these numbers to 9 (with 
massless leptons) and 12 including masses. The conclusion is that there is a definite 
number of independent experimental observables that can be extracted from the angular 
analysis, and this number riexp is known from symmetry arguments. From a pragmatic 
standpoint, the symmetry formalism can be substituted by this set of riexp independent 
observableqj- This set of riexp independent observables can be considered complete, in the 
sense that any additional observable can be expressed as a function of the observables in 
this set. Such a complete non-redundant set may be conveniently called a basis. 

In this context, the natural question is what is the best choice for the observables in 
the basis. The answer is not different from before: these observables should satisfy the 
desired properties of reduced hadronic uncertainty and good sensitivity to New Physics. 
The observables Ap^'*'^\ v4f''''°\ h^^'^'^'^'^) proposed in Refs. [H El [3 El [9] are excel- 
lent candidates, since they were designed to satisfy these requirements. The question is 
whether this set of observables contains a basis, or if other observables must be introduced. 

The purpose of this paper is to give an answer to this issue and to construct, in a 
systematic way, an optimal basis of observables related to the angular distribution of 
the decay — )■ K*^{^ K~^7T~)i~^i~ , including lepton masses and scalar contributions 
(and valid in the presence of tensor operators). We do not consider in this paper the 
CP-conjugated mode nor the corresponding CP asymmetries. Anticipating some of the 
results of the paper, we shall see that: 

■^Evcn if the term 'symmetry" usually denotes an invariance of the Lagrangian, in the present paper 
and following Ref. [7] , it refers to an invariance of the angular distribution under a rotation in the space 
of spin amplitudes. 

Still, we will use the symmetries at some points to furnish more formal proofs of certain aspects of 
the approach. Also, the symmetry formalism can be used, for instance, to obtain explicit expressions of 
the transversity amplitudes in terms of the coefficients Ji. 
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1. The optimal basis contains two types of observables: observables with LO depen- 
dence on the soft form factors (Form Factor Dependent (FFD) observables) and 
observables free from this dependence at LO (Form Factor Independent (FFI) ob- 
servables). The FFD observables suffer from large hadronic uncertainties due to the 
poorly known soft form factors. For this reason, the goal is to maximize the number 
of FFI observables in our basis. Examples of FFD observables are Ji{q'^), or Aim. 
We will choose as FFD observables for our basis the differential rate dV /dq^ and the 
forward-backward asymmetry ApB, although other choices are also possible. 

2. In the case of massless leptons, the optimal basis according to the previous counting, 
must contain 8 independent observables: two FFD observables {dV /dq^ and A^b) 
and six clean FFI observables. Examples of FFI observables are Auj^ ^ , Auj^ 

or if^^'^'^'^'^-*. FFI observables can be constructed in a systematic way as we will 
show in Section [3l leading to a set of observables -Pi,2,3,4,5,6; which we call primary 
observables. Pi, 2,3,4,5 are directly related to all already known observables, but is 
new, and it is necessary for obtaining full information from the angular distribution. 

3. In the massive lepton case, the basis must contain 10 observables (2 FFD and 8 
FFI). Eight of those observables are dT/dq"^, Ap^ and Pi with i=1...6. The two 
remaining FFI observables, which vanish in the massless limit, will be called Mi 
and M2 and have never been considered before. 

4. In the scalar case with massive leptons the counting of symmetries establishes the 
existence of 12 independent observables. The full set of observables is composed by 
the two FFD dV/dq'^ and Afb together with 10 FFI Pi,2,3,4,5,6, Mi, M2, ^1 and ^2. 
The two new observables 5*1 and S2 vanish in the absence of scalar contributions. 

5. Any conceivable FFD observable can be written as a function of dT/dq"^, Ap^, Pii^t 
and Si, but most importantly, any conceivable theoretically clean (FFI) observable 
can be written as a function of the Pj, Mi and Si alone. At the same time, each of 
these observables contains unique information. It is in this sense that the basis is 
optimal and complete. 

The structure of the paper is the following. In Section[2]we briefly review the symmetry 
formalism for the angular distribution of 3^ — )■ A'*°(— )■ K7i)£~^i~. In Section |3] we build 
the basis of observables in the massless case, and show how the full set of observables that 
have been considered in the literature can be recovered from this basis. We shall keep mass 
effects at this stage to make the generalization in Section |4] most straightforward. After 
presenting the generalization to the massive case, and introducing the massive observables, 
we write the full angular distribution in terms of the observables in the basis. 

In Section [5] we include the effect of scalar operators, and show how the previous 
results are modified. In particular, we introduce two extra observables, 5*1 and S'2, that 
vanish in the absence of scalar contributions. In Section |6] we present the most general 
expressions for the three uniangular distributions in terms of the observables in the basis. 
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In Section[7]we study the New Physics sensitivity of the proposed observables. For that 
purpose we study the SM contribution including NLO effects using QCD Factorization, 
hadronic uncertainties and an estimate of A/m?, corrections. We also consider how the 
SM expectations are modified in several NP scenarios. We analyze the position of the 
zeroes of the observables as well as those NP scenarios that affect most strongly each of 
the observables. 

In Section [8] we summarize the relevant results of the paper. Finally, the core of the 
mathematical machinery related to the symmetry formalism, including constructive proofs 
of existence of the continuous symmetries has been collected in Appendix[Xl This appendix 
contains the explicit form of the symmetry transformations among the amplitudes in the 
massless (appendix lA.ip . the massive (appendix IA.2P and scalar case (appendix I A. 30 . In 
appendix [B] we present the building blocks of the observables in the large recoil limit. 



2 Symmetries of the angular distribution 

The coefficients of the distribution given in Eq. ([1]) can be written in terms of transversity 
amplitudes. In the massless case there are six such complex amplitudes: A^'^, A^'^ and 
A^'^. An additional complex amplitude At is required in the massive case, and in the 
presence of scalar contributions a new amplitude As must be included. The expressions 
for these coefficients read, 

j^^ = (1±M. + |Af |2 + \Af\^ + \Af\^] + ^Re (AiAf + A^Af) , 



Jl, = \A^\' + + ^\\At\-' + 2Re{A^A^^)\ + \A 



Ami 



L aR* 



J^^ = ^ [|Ai|2 + |Af|2 + |A^|2 + |Af|2] 



1 



^2c = -/3,Mi^or+i<r] 
1 



Js = V2Pe he{A^Ai* - A^A^) - ^ Re(A[ A* + AfAs 



Jes = [Re{AlAi* - AfAT)] , 



J,, = 4Pe^MAf;A*s + A^*As) 



Jl = V2/3£ 



L aL* aR aR*\ I 



lm{Al;Af - A^A 



lm{AiA*s - AT As)) 



Js = -j=f3! [lm(AoMi* + A^A^ 



Jg = /3f [Im(4*4 + AfAf)] , (3) 
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where the parameter is given by 




/3. = Wi--:/. (4) 



We will distinguish between the three cases of interest: massless leptons and no scalar 
amplitude, massive leptons and no scalar amplitude, and massive leptons plus a scalar 
amplitude. We also show that tensor contributions cannot change the picture. 

A. Massless leptons, no scalars 

In this case we can put As — )■ 0, drop the mj terms and set /?£ — )■ 1. There are in total six 
complex transversity amplitudes: A^'^, A^'^ and A^'^, which add up to 2nA = 12 real 
theoretical quantities. However, an infinitesimal transformation of the distribution (see 
[7]) shows that there must be = 4 continuous transformations between the A^ ' that 
leave invariant the angular distribution. Two of them are simple phase transformations: 
Af — > e^'^^Af and Af — >■ e*'^^Af, while the other two mix L and R amplitudes (see 
Appendix [A] and Ref. [7] for the explicit form of these transformations). According to 
Eq. (I2l), there must be precisely riexp = 8 independent observables. This implies, in turn, 
that there should be 4 relationships among the 12 coefficients Ji{q'^)- Three of them are 
straightforward: Jgc = 0, Ju = 3J2s and Jic = —J2c, while the remaining one is more 
involved [7j: 

^ _ ^ (2 J2. + J3) (4J| + Pp^) + {2J2S - J3) {PjJl + 4J|) 

J2c — — ^ 



16J|,-(4J| + /32J|, + 4J|) 

16 J|. - (4 J| + PIJI + 4 J|) ' ^ ^ 

where Pe 1 should be understood in the massless case. The derivation of this expression 
requires the symmetry formalism and will be outlined later. 



B. Massive leptons, no scalars 

In this case we just set — > in Eq. ([3]). Now there are seven complex transversity 
amplitudes, including At, which add up to 2n^ = 14 real theoretical quantities. As dis- 
cussed in Appendix IA.21 there are = 4 continuous transformations that leave the J, 
invariant (these are different from the symmetries of the massless case, since those are 
broken by mass effects). Two of these symmetries are phase rotations: At — ?■ e^^^At and 
^o'li ± ~^ ^*'^^o '|| ±' "^hile the other two are nonlinear transformations (see Appendix lA. 21) . 
According to Eq. (|2]), there must be precisely rigxp = 10 independent observables, which 
means that 2 relationships between the coefficients Jj can be found. The relationships 
Jis = 3J2S and Jic = —J2C are no longer satisfied; however, Jgc = and Eq. ([5]) remain 
exactly true in the massive case (this was the reason for keeping the factors /?£ explicit in 
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Eq. ([5])). Notice that in Ref. [7] the discussion was hmited to the massless lepton case, so 
Eq. ([5]) generahzes the relation in Ref. [7] to the massive case. 

C. Massive leptons plus scalars 

In this case we deal with 8 complex transversity amplitudes, which add up to 2n^ = 16 
real theoretical variables. As demonstrated in Appendix IA.3I (see also Ref. [7]), there are 
Us = 4: symmetries. This means there must be exactly rzexp = 12 independent observables, 
which implies that all the J, are independent, and none of the previous relations hold in 
this case. 

D. Massive leptons, scalars and tensors 

The fact that all the Jj are independent in the scalar case can be used to go a bit further 
in the reasoning. Imagine we wanted to include NP contributions from tensor operators. 
Then we would expect, at least, a new amplitude At modifying somehow the angular dis- 
tribution [Eqs. ([3])]. However, according to Eq. ([2]) and since riexp is as large as it can be, 
for each new amplitude At there must be two extra symmetries. These symmetries must 
disappear in the hmit At — 0, and therefore they can be used to set At — >■ 0. Another 
way to see this is the following: since the Jj are all independent in the scalar case, one 
can always obtain a set of A^'^'q, A[, A'g that reproduce the angular distribution in the 
presence of new amplitudes such as At- Therefore, new tensor operators can only give 
new contributions to existing amplitudes, meaning that the basis of observables defined 
in the scalar case remains unchanged in the presence of tensors. 

We will now discuss in turn the relevant set of observables that one can consider in 
each of these three cases of interest. 



3 Observables for massless leptons 

Not any observable constructed from the transversity amplitudes can be obtained from the 
angular distribution. As a necessary and sufficient condition, such an observable must be 
invariant under the symmetry transformations of the transversity amplitudes A's; we then 
say that the observable respects the symmetries of the angular distribution. Fortunately, 
there exists a systematic procedure to construct all such possible observables. 
We start defining the following complex vectors [7], 
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With these vectors we can construct the products |njp 

R\2 _ ^"^25 "~ J3 



n\ni and n\ rij, 



\nii 



PI 

2J2S + J3 



6s 



2iJa 



nl n 



Pf. 



n 



t „ , _ al* al _ ar ar* 



2/3| 
2 J4 — ifi^J-j 

J5 - 2zJ8 



(7) 



These quantities automatically respect the symmetries of the angular distribution, since 
they can be expressed in terms of the Jj. Considering real and imaginary parts, there 
are 9 real quantities that encode all the information of the angular distribution, and by 
combining them one can construct systematically all possible allowed observables consis- 
tent with the symmetry requirements. However they are not all independent: any set of 
complex 2- vectors {no, ny, nj_} satisfies 

|(njn_L)|noP - (nj no)(nJ, n^) |^ = (|non?^||P - \nln\\\^){\no\^\nxf - \nlnif) . (8) 

Using Eqs. ([7]), this relation translates precisely into the relation for the Jj given in Eq. (|5]). 

Now that the formalism assures the systematic construction of observables that respect 
the symmetries of the angular distribution, we must focus on the cancellation of hadronic 
form factors. At leading order in l/mh and a^, and at large recoil {Ek* — )■ 00), the 
transversity amplitudes Aq'^, A^'^ and A^^ can be written as: 



A^'^ = V2NmBil - s) 



WJ 



A 



L,R 



A 



L,R 



■V2NmBil - s 
NruBil - sf 



c 



K* 



(Cf - Cf) T (Cio - C[,) + 2rn,{Cf - C 



cff 



K" 



(9) 



where s = Q'^/m^, rfij = rrii/mB, and terms of O^rhj^,) have been neglected. The normal- 
ization is given by 



3 ■ 2i%5^| 



(10) 



with A = [q"^ — (m^ +mj^.)^][g^ — (m^ — mx*)^]- Therefore, at first order, we have no oc 
and ny , n_L oc This establishes a clear guideline in the construction of clean observables, 
as ratios of quantities in Eq. where the C,\\,± cancel [Form Factor Independent (FFI) 
observables] . 

Before providing a complete list of observables constructed according to this procedure, 
we should note the following. There are 8 independent quantities in Eq. (JTj) that constitute 
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the building blocks of the observables. The soft form factors ^\\^± can be thought of as 2 
irreducible normalization factors in the products nlrij, and therefore one cannot construct 
8 independent combinations where the soft form factors cancel. The best one can do is 
to construct 6 clean observables, plus 2 observables that contain the information on the 
two form factors — or Form Factor Dependent (FFD) observables-. 

For these two FFD observables we can choose, quite naturally, the angular-integrated 
differential decay rate dT/dq"^, and the forward-backward asymmetry ApB- 

— = / dcoseidcoseKd(f)-—— ^—z ^-7- = 7 (3 Jic + 6 Ji, - Jac - 2 , (11) 

dq^ J dq^ rfcos Ok "Cos Oi dcp 4 



A 



FB 



dV/dq^ 



dq^dcosOi 3Jic + QJu - J2C - 2J2s 



Notice that, while Eq. flTT]) is completely general, in the last equality of Eq. (1121) we have 
assumed that Jgc = due to the absence of scalar contributions. In the massless case, 
since Ju = "iJ-zs and Jic = — J2C5 these expressions simplify to dV /dq^ = Jic + 4J2s and 

Afb = -3J6./[4(Jie + 4J2.)]. 

For the six (clean) FFI observables we choose the following set: 



\n± 


2 _ 


nil 


2 T 


\n± 


2 + 


n|| 


2 2J2. 



^1 = = ^ ^ (13) 



Re(n^i nil) „ Jqs 



_ Im(n^^nii) _ Jg 

- |n||P + |nJ2 - -4J2: ' ^^^^ 

p, = 4£hi = ^^^=, (16) 





-\/"~^2c(2J2s ~ 


-^3) ' 


Re(no n±) 








1^— 2 J2c(2J2s 


+ J3) ' 


Im(nJ n||) 







although other similar ratios are possible. We have used the following criteria for choosing 
among the different possible FFI observables: (1) they are simple ratios of the quantities 
in Eq. ([7]) where the form factors ^±^\\ cancel, (2) they take values in the range [-1,1], and 
(3) they show good sensitivity to selected New Physics (see Section [7j). To summarize, 
the complete basis of observables in the massless case is given by: 

Omi=o = {^^' ^FB, -Pi, P2, P3, Pi-, P5, -Pej (19) 
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where the six Pi are clean observables. 

As mentioned, all possible observables can be expressed in terms of the observables in 
the basis. In particular, all known observables can be related to this set. For example, 
the usual FFI observables can be expressed in terms of -Pi,2,3,4,5,6: 



.(2) 



Pi , 




- P2 - 4P| - 


4P2 




A?-") = -2P3 






Pa. 








(5)n2 


h + (Ap^)' 


+ = 


1, presented in 



Ap^ = 2P2 , 4"") = -2P3 , (20) 

Also, the relationship {2A\ 
follows trivially in terms of the Pj. 

In the case of the observables A^'* and A^\ the corresponding expressions are more 
involved. They can be expressed in terms of our basis by first recovering their expression 
in terms of the J's: 



^ ^/ -2J2c(2J2. + Ja) ' ^ yUi + Ppy ^ ' 



and substituting the J's in terms of our basis of observables as given in Eqs. ( l3T]) -( H0l) in 
the following section. 

The same can be done for all the known FFD observables, which can be expressed in 
terms of Pj, dV /dq^ and Afb, for example: 

P^ = l-P, = -^^ (22) 
^ ^ 3 P2 ^ ' 

Aim = - (23) 

4 Observables in the massive case 

All the coefficients in the massive case can be expressed in terms of the quantities of 
Eq. ([7]), with the exception of Jic and Jis, which can be written as: 

^ fi 1 2 I 1 2~l ^TTh f r 'T' \ ^ T' t=(=~l/\ 

■Jis = — ^ — [\n±\ + \n\\\ \ + |_n|| criny + n^^^aiu^^ - nj_ain± - n'j_ainj_\ , (24) 

Jic = \no\^ + —if [2| AT + n^(^ino + njcnn*] . (25) 
where cxi is the Pauli matrix: 

The important point demonstrated in Appendix IA.2I is that in this case the symmetries, 
like in the massless case, can be expressed as a single unitary rotation U on no,±,||. This 



10 



means that all the products [Eq.Q] that were invariant in the massless case are still 
invariant under the new symmetries (so all the observables defined in Section |3] are still 
valid), which in turn means that the new terms appearing in Jic and Ji^: 

^Re [AlAf + A\Af) , ^[\A,\' + 2Re(A^<*)] , 

must be invariant by themselves. This is a key idea in order to find the continuous 
symmetries (see Appendix IA.2|) . but it is also crucial to the construction of the new 



observables, the obvious ones being: 

T t 

2m? 1 niairiw^nlain 



2mj 




t * 


PjJu-{2 + Pj)J2s 




fif |n|||2 + |nx|2 








2 p + riQairiQ + n\ainl 


I^JJlc + J2C 













M2 = ^ ■ ' " ' V 9 = -^^^^^^^-^^ . (2 



We note that these observables are of the FFI type, and thus theoretically clean. This 
can be inferred from the large recoil limit expressions ([9]) and 

A, = ^^^(1 - sf [2(Cio - C[,) + ^{Cp - C'p)] i\\{EK^) . (29) 



Moreover, Mi and M2 vanish in the massless limit (from the right hand side of Eqs. (|2j 
and fl28|) one can see that this follows from the relationships of the massless case, Jis = 3J2s 
and Jic = —J2C1 that are broken for non-zero lepton masses; in fact Mi measures the 
breaking of the relation Jis = 3J2s, while M2 measures the breaking of Jic = —J2c)- 

From the discussion in Section [21 together with the observations that Mi and M2 
vanish for — 0, and that the observables of the massless case are still valid, one 
concludes that these two observables complete the basis of 10 independent observables of 
the massive case. This basis is: 



Vn.^o = ^FB, Pi, P2, P3, Pa, P5, Pg, Ml, M2} . (30) 



The coefficients Jj of the angular distribution are themselves observables (of the FFD 
type), and it is interesting to express them in terms of dT/dq'^, ApB and Pi. The explicit 
expressions read: 



2 + (l + 4Mi)/3| dV 

= We — ' 



= ...TAfL^, (6/5^ + [3 + (1 + 6Mi);3,^] x) ^ ' (^2) 



3A(3 + 3M2 + /3|) ^ ■ L ■ V . ^^^^^^^(ig2 

^ dT 
6 dq'^ 



J2s = -—X-To , (33) 
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= -^ITTMkTW) ^ " ^ 1^ . (34) 

1 cir 

■^3 = ^Xt^ , (35) 




^4 = ^^/ oVo./^,^2 (6/3^ + [3 + (1 + 6Mi)/3|] X) X 1^ , (36) 



^5 = ^J o^, 3M, + P- ^^^' + [3 + (1 + 6MO/3|] X) X |i , (37) 



4P2 , , 

= , (38) 



3 dq 




J. - J^^^^ (6/5^ + [3 + (1 + 6M,)/3|] X) X |; , (39) 



where, in this case, 

X = ■ (41) 

-'2 

The observable x is well defined because ApB and P2 share the same zeroes [see Eqs. (1121) 
and f|T^ ]. This will change after including scalars, but also x will change. On the other 
hand, since Ji^ > J2C [as can be checked from Eq. Q], it follows that 3 + 3M2 + /3| > 0, 
and no vanishing denominators can occur in Eqs. ( !3T|) - (l40|) . 

Notice that Jg is absent from this list, as it is not an independent coefficient in the 
absence of scalar contributions. In order to obtain its expression in terms of the observ- 
ables one should write Jg as a function of the other Jj using Eq. ([5]), and then plug in 
the expressions ( 1311) -( 140|) . Alternatively, one might want to leave Jg as a free parameter 
when fitting the angular distribution in terms of the observables, and then check the 
relationship in Eq. to look for scalar contributions (see next section). 

Using the substitutions Mi,M2 — ^ and — )■ 1, Eqs. (13T1) - P0l) transform into the 
corresponding massless case expressions. 



5 Inclusion of Scalar Operators 

In the presence of scalar operators, a new amplitude As appears. This makes 16 the 
number of real theoretical degrees of freedom (8 complex amplitudes): 2n^ = 16. There 
are 4 symmetries, making rzexp = 12. This means that all the Jj are independent, which 
in turn implies that the relationship Jqc = and Eq. ([5]) do not hold any longer. This 
observation allows us to choose the two extra observables that are needed in the basis 
in the presence of scalar contributions, as the amount by which these relationships are 
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broken. This provides two independent observables, Si and 5*2, that vanish in the absence 
of scalars. The first observable is: 

Si = , (42) 

Ame J2c 

which measures the breaking of the relation Jgc = 0. The second one is: 

^2 



lQJ2cJl - 4^2cJ| + lQJ2sJl + SJsJi + 16J2.^| - SJsJi + 1674^8-^9 " 4J2c^| 



/^scalars 



J 2c J 2s 

n2^'^2sJ'i ~ "^JsJi ~ 4J4J5J6S ~ J2cJqs + 4J2SJ7 + ^J^J^ — AJqsJ^Js — AJ^J-fJg , , 
+ P J ^2 5 \^'^) 

J 2c J 2s 

which gives a measure of the violation of Eq. (|5]). One can easily check that both observ- 
ables are of the FFI type, by noting that the large recoil expression for As is: 

As = -2 ^(1 - sf [Cs - C's] ^iiiEK,) . (44) 
rriK* 

While most of the results in the previous sections remain unchanged in the presence of 
scalars, some differences must be clarified. When As 7^ 0, the observables P5, Pq and 
M2 get modified. In particular Eqs. (fT7|) . ( !T8|) and f l28|) in terms of the vectors rii do not 
hold since new terms proportional to As arise, and these observables must be redefined. 
The simplest way to generalize these observables in presence of scalars is simply to use 
their definition in terms of the Ji in Eqs. ( IT7|) . (fT8|) and (!28|) . These three observables 
are, together with 5*1 and S2, the only ones sensitive to scalar contributions. With this in 
mind, an optimal basis of observables in the presence of scalars reads: 

dT 

In the case of ApB, we keep the definition in Eq. (fT2|) as the angular integral, which means 
that now: 

^ 3J6g + 3J6c/2 ^^g^ 

3Jic + QJis — J2C ~ 2J2S 
Because of this, Eqs. (!3T]) -( l40|) are modified slightly. First, when S'j 7^ we have: 

(3 + 3M2 + /3|)Afb + (m,/v^) Q^eSi , 

K = 1= (47) 

(3 + 3M2 + /3|)P2 - (m,/v^)(3 + 6/3|Mi + /32)5i 

in (13T])- (|40|) . Furthermore, Jec is non-zero: 

8m,5i (6/32M1 + /3| + 3)Afb + 6/3,P2 

Jf\r = 1=^ ;= TTT • (4o) 

3 (3M2 + /3| + 3)P2 - (m,/ v^) (6/3|Mi + /3| + 3) ^1 dq' 
The generalization of x iii Eq. (1471) could be expected: in the case — t- 0, the zeroes of 
Ayb and P2 coincide, because they are both proportional to Jqs [see Eqs. f|T2l) and f lT^ ]. 
and that is the reason for which the combination ylFB/P2 is well defined in Eqs. ( l22l) . 
( 123|) and (13T]) -( 140|) . However, turning on As moves the zero in A^b away from that of P2, 
making ylFB/P2 singular when Jgs = 0. On the other hand, the expression in Eq. (147|) is 
regular at every point and goes to x = ^fb/P2 only in the limit As — )■ 0. 



|^,AfB,Pi,P2,P3,P4,P5,P6,Mi,M2,5i,52} . (45) 
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6 Uniangular distributions 



A full angular analysis designed to extract the complete set of g^-dependent observables 
requires a good deal of statistics, and will be possible at LHCb not before an integrated 
luminosity > 10 fb~^ has been collected. However, certain angular observables are avail- 
able from partially integrated distributions, and experimental analyses have focused to 
this day on uniangular distributions, leading to the set of measured observables reviewed 
in Section [H 

Starting from the full angular distribution of Eq. ([1]), the three uniangular distributions 
can be obtained: 

d'^T 1 

' (6 Jis + 3 Jic — 2 — + 4 J3 COS 2(f) + 4 Jg sin ! 



d^T 
dq^ dcos 9e 

d^T 



Svr 
1 



QJu + SJic + (6J2S + 3J2c) cos26'^ + (6J6s + 3J6c) cosOi 



{9Ju - 3 J2.) sin^ 9k + (9 Ji^ - 3 J2c) cos^ Ok 



(49) 



(50) 



(51) 



dq"^ dcos 6k 

Substituting the expressions fl3T|l - fj4Ql) for the Jj coefficients, the uniangular distributions 
in the presence of scalars can be written as functions of the observables as follows: 



dq^dcl 



1 

2^ 



1 - y cos 20 + sin 20 



dq^ 



(52) 



dq^ dcos 9 1 



3(M2 + 1) 



-A 



2(3M2 + /3| + 3) 
2(3M2 + (5j + 3) 



+ 



4/3|Mi + M2 + /3| + 3 
8(3M2 + (3j + 3) 

3/3£(4/3|Mi + M2 + /3i 



X - v4fb cos On 
' + 3) 



(53) 



8(3M2 + /3| + 3) 



X cos29i 



dr 

dq^ 



d'^V 



dq^ dcos 9k 



3 6/3|Mi + /3| + 3 \ . 6/3|Mi + /3| + 3 . 



dT 
dq^ 



(54) 



where x is given by: 

' ^FB 



X 



if A, 



(3 + 3M2 + /3|)Afb + (m./y/^) 6/3,^1 
I (3 + 3M2 + /3|)P2 - (m,/v^)(3 + 6/3|Mi + Pj)Si 



(55) 



if A^T^O 



Note that in the limit — all scalar effects disappear. The corresponding well-known 
expressions for the uniangular distributions in the massless case are obtained by setting 
0, /3f ^ 1 and Mi, M2 0. 
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C^ifJ-b) 


CaifJ^b) 


CeifJ-b) 


Cfif^b) 


Cfif^b) 


CgifJ-b) 


Cio(/"fe) 


-0.2632 


1.0111 


-0.0055 


-0.0806 


0.0004 


0.0009 


-0.2923 


-0.1663 


4.0749 


-4.3085 



Table 1: NNLO Wilson coefficients in the Standard Model at the scale fib = 4.8 GeV [27]. 
In the computation of the observables, we consider a variation of /x G [/ib/2, 2/1^]. The 
coefficients Cg and Cg^ are related though Cf^ = Cg + F(g^) (see Ref. [29]). 



7 New Physics Sensitivity of the Observables 

In this section we analyze and discuss the New Physics sensitivity of the full set of ob- 
servables 0^^}^Q^. In particular, we study the impact of New Physics contributions to the 
Wilson Coefficients: 

C, = Cf *^ + 5Ci , (56) 

where i = 7^'\9^'\l0^'\ S, P, always taking into account the existing bounds from other 
processes that constrain substantially the New Physics contributions 6Ci. 

We consider the 10 FFI observables -Pi,2,3,4,5,6; ^1,2 and 5*1^2 in terms of the transversity 
amplitudes A^^'^, A^^, At and As- These amplitudes can be written in terms of the 
Wilson coefficients Ci and a set of seven form factors V{q'^), /lo,i,2(9^) and Ti^2,3(9^) (see 
for example Ref. [21]). We ffist consider the SM contribution to the observables including 
NLO corrections, hadronic uncertainties and an estimate of the A/rub corrections. We 
will see that indeed these FFI observables show reduced hadronic uncertainties. 

After having the SM contribution under control, we consider NP contributions in sev- 
eral different scenarios, all of them compatible with current bounds from other processes, 
and study the possible deviations from the SM. The outcome of this analysis is shown in 
Figs. [Ha and Tables EP 

7.1 SM contribution and hadronic uncertainties 

The SM Wilson coefficients at the matching scale fiQ = 2Mw, and their running from fiQ 
down to /ifo = 4.8 GeV, as well as the running of quark masses and couplings, are obtained 
following Refs. [221 |23l EH [251 [26] (see also Ref. [27]). For reference we quote in Table [D 
the used values for the Wilson coefficients at the hadronic scale, taken from Ref. [27] . 

Concerning the seven B — )■ K* form factors (y(g^), ^0,1,2(9^) and Ti_2,3(q'^)), their 
g^-dependence is parametrized following Ref. [28], giving more conservative uncertainties 
than other parameterizations. Their values at = are given in the same reference, 
obtained from light-cone sum rules with B distribution amplitudes. The definitions of 
the soft form factors C,\\,± in terms of the full form factors are given in Ref. [30]. All 
numerical inputs used in this analysis are the same as the ones tabulated in Section 2.1 
of Ref. [27]. Maximum and minimum values for ^||(g^) and ^_L(g^) give rise to the grey 
regions in Figs. [T|2] around the central SM value. 

At this point, the rest of the hadronic uncertainties are calculated: 
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1. The renormalization scale /i6 is varied between 2.4 and 9.6 GeV. 

2. The value of rhc = mc/m^ is varied in the range rhc = 0.29±0.02, according to the 
discussion in Refs. [291 131] . 

3. The uncertainty related to the factor that determines the relative size of the hard 
scattering term vs. the form factor, (see Eq. (55) of Ref. [29] and discussion below), 
is estimated at the level of a 30%. 

These uncertainties are added in quadrature together with the uncertainties related to 
the form factors, giving rise to the orange bands in Figs. [T][2| on top of the gray bands 
(which include only the form factor uncertainties). 

As a third step, A/rrih contributions are estimated following the procedure in Section 
2.3 of Ref. [7], but widening the error band to include a 68.2% of the probability (as 
opposed to the 66% used in that reference). This uncertainty is added in quadrature to 
the rest of the uncertainties computed before, giving rise to the light green bands (5% 
A/rrib correction) and the dark green bands (10% A/rrib correction) in Figs. [T1I21 

7.2 New Physics 

The impact of NP on -Pi, 2,3,4,5,6; ^1,2 and 5*1,2 is shown in Figs. [T|2l According to the 
model-independent fit of Ref. [27] (updated in Ref. [32]), three sets of values for Cj are 
chosen in order to represent the NP impact on the observables. 

Scenario A {Ct^C'j): In this scenario, Cj and are chosen according to the allowed re- 
gions obtained in the analysis of Refs. [271 132], while the rest of the coefficients are set to 
their SM values. In particular, Cj is set to values where deviations between experimental 
data and the SM are maximal. We choose two subscenarios: Scenario A.l corresponds to 
a "SM-like" point in the Cj — C'-j plane belonging to the allowed connected region where 
the SM lives (see Fig. 2 of Ref. [32]). Scenario A. 2 corresponds to a "non SM-like" point 
belonging to a disconnected region, most likely to be probed in the near future by im- 
proved measurements of the branching ratio of i? — )■ Xsfi~^fi~ ■ The values of the relevant 
Wilson coefficients in these scenarios are summarized in the first column of Table [21 

Scenario B {Cg,Cio): In this scenario, C7 and Cj are fixed with small NP contributions, 
while Cg and Ciq take maximum allowed values compatible with the chosen C7 and Cy. 
We distinguish between Scenario B.l, where only Cg receives a non-zero NP contribution, 
and Scenario B.2 where the NP enters only in Ciq. The values of the relevant Wilson 
coefficients in these scenarios are summarized in the second column of Table [21 

Scenario C {Cq,C[q): In this scenario, C7, C'j are fixed as in Scenario B, and Cg, Ciq are 
SM, whereas Cg and C[q take the maximum allowed values compatible with the given C7, 
C7, according to Refs. [271 [32]. We again distinguish between Scenario C.l, where only Cg 
receives a non- zero NP contribution, and Scenario C.2 where the NP affects only to C[q. 
The values of the relevant Wilson coefficients in these scenarios are summarized in the 
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Scenario A Scenario B Scenario C 





A.l 


A.2 


B.l 


B.2 


C.l 


C.2 




-0.041 


0.25 


-0.002 


-0.002 


-0.002 


-0.002 




-0.114 


-0.414 


-0.006 


-0.006 


-0.006 


-0.006 


6Cg{j2b) 






-1.25 


















-3 












3 


















-3.5 



Table 2: Wilson Coefficients at the hadronic scale /Xb = 4.8 GeV within Scenarios A,B,C. 





Scenario D.l 


Scenario D.2 


Scenario D.3 




0.1 + 0.5i 


1.5 + 0.3Z 




SCg{fib) 


-1.4 


-8 + 2i 




5Cio(/i6) 


1 - 1.5i 


8-2i 




SC'^{fib) 






-0.3 -0.1 z 








3 + i 








-0.6 + 2z 



Table 3: Wilson Coefficients at the hadronic scale fib = 4.8 GeV within Scenario D. 

third column of Table [2l 

Observables P3 and Pg are mostly sensitive to imaginary components of the Wilson co- 
efficients, since they are built out of imaginary parts of amplitude products [see Eqs. ( fT5i) 
and (ITSl) ]. In order to test this dependence, we consider a fourth scenario with complex 
NP contributions to the Wilson coefficients Cy^g^io and C'r^gio- 

Scenario D (complex WC's): In this scenario, the NP contributions SC^'^ 10 take com- 
plex values. We consider three possibilities. In Scenarios D.l and D.2, only Cy^g^io receive 
NP contributions. Scenario D.l consists on a point inside the "SM-like" allowed region 
found in Ref. [33], while Scenario D.2 is a point in the other "non SM-like" region of 
Ref. [33]. In Scenario D.3, only Cjg iQ are affected by NP. The values chosen for the Wil- 
son coefficients in these scenarios are summarized in Table [31 

Finally, we consider two additional scenarios with scalar and pseudoscalar New Physics 
contributions, to study the scalar observables 6*1^2 and the pseudoscalar sensitivity of M2: 
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Scenario S Scenario P 





S.l 


S.2 


S.3 


S.4 


P.l 


P. 2 


P.3 


P.4 


{Cs - C's){fit) 


0.03 


0.01 


-0.01 


-0.03 










{Cp - C'p){fi,) 










0.07 


0.0467 


0.0233 






Table 4: Wilson Coefficients at the scale /ij, = 4.8 GeV within Scenarios S and P. 



Scenario S (Cs—C'g): All Wilson coefficients are SM except for Cg\ Since the amplitudes 
are only sensitive to the difference Cs — Cg, we consider four different values for this 
difference, all compatible with the latest Bg — )■ yU"*"//" bounds [M] : 

BR{B, < 4.5 X 10"^ (57) 

at 95% confidence level. These four values constitute Scenarios S.l to S.4 and are sum- 
marized in Table |H 

Scenario P (Cp — C'p): In this case, all Wilson coefficients are SM except for Cj^^ A gam, 
since the amplitudes are only sensitive to the difference Cp — Cp, we consider four different 
values compatible with the Bg — )■ ^~ bound. These four values constitute Scenarios 
P.l to P.4 and are summarized in Table |H 

The set of observables can be divided in two groups: -Pi,2,3.4,5,6 and Mi, which are 
only sensitive to Cj with i = 7^'\ 9^'\ 10'-'^ constitute the ffist group. In the second group 
we include M2 and Si 2, which are also sensitive to Cs and Cp. In principle, P5 and Pg 
contain also Cs [see discussion below Eq. flH|) ] but the overall sensitivity, considering the 
present bounds on Cs, is negligible (for this reason we do not present the curves for P5 g 
in Scenario S). We will focus on the case £ = /i in all considerations concerning lepton 
mass effects. 

Within the ffist group we have the observables Pi, P3 and Pq, that are suppressed in 
the SM in all the region, but which can take sizeable non-vanishing values in specific 
NP scenarios, as shown in the left column of Fig. [1] 

The rest of observables in the ffist group, namely P2, P4, P5 and Mi are no n- vanishing 
aheady in the SM and present a non-trivial g^-dependence. They all contain a zero 
at a value of within the experimentally accessible region 1-6 GeV^. At LO, these 
zeroes occur at the positions specified in the left column of Table [5] (these include NLO 
corrections in the Wilson coefficients). At NLO, hadronic corrections shift the zeroes by 
amounts that can be computed in QCD factorization (see the second column in Table [5l). 
In the presence of NP, the positions of the zeroes are substantially modified (moved to 
lower or higher values) as can be seen in Figs. [T1I21 In the third column of Table |5l we 
summarize the position of the zeroes in the most relevant NP scenarios considered above. 
In some NP scenarios, the zero can even disappear from the low-g^ region. This is the 
case for example for P4 in Scenario A. 2, or P5 in Scenario C.2. 
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Figure 1: Observables -Pi, 2,3,4,5,6 in the SM at NLO, including all hadronic uncertainties 
(wide bands) as explained in Section 17.11 The solid and dashed curves correspond to 
the NP scenarios exposed in the text: Scenario A.l (blue solid), Scenario A. 2 (blue 
dashed), Scenario B.l (red sohd). Scenario B.2 (red dashed). Scenario C.l (brown solid) 
and Scenario C.2 (brown dashed). Scenarios D.l, D.2 and D.3 are explicitly indicated in 
P3 and Pq. The Wilson coefficients responsible for the largest deviations are highlighted. 
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Figure 2: Observables Mi^2 (for i = fi), dT/dq^ and ApB- The bands correspond to 
the SM at NLO and with all hadronic uncertainties. Notice that Si and S2 are strictly 
zero in the SM. Dashed and solid lines in Mi, dV /dq^ and Ayb correspond to NP Scenarios 
A,B,C as in Fig. [H In the case of M2 and 51^2, the curves correspond to: Scenario S.l 
(blue solid), Scenario S.2 (blue dashed). Scenario S.3 (red dashed). Scenario S.4 (red solid). 
Scenario P.l (gray solid). Scenario P.2 (gray dashed). Scenario P.3 (magenta dashed) and 
Scenario S.4 (magenta solid). The Wilson coefficients responsible for the largest deviations 
are indicated. 
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Obs. 


Qq'^'^^ at Large Recoil 


g2SA/ (NLO) 




P2 


2mi,m B cf 


3.06 


3.93 


5.23 [B.l] 


9 


P4 


2mbmBimB Cf + 2mbCf) Cf 


1.58 


1.87 


3 75 \C. 21 
3.55 [B.2] 


2m, + m D i'r'efr2 I n2 \ 


P5 


mbm\ ^7^ 


1.64 


2.23 


3.25 [B.l] 
3.03 [A.I] 


rUbCf + mB Cf 


Ml 


2mfem B cf 

L-g — C-IO 


1.61 


2.07 


3.15 [B.2] 



Table 5: Position of the zeroes of observables -P2,4,5 and Mi, in the SM at large recoil, 
at NLO, and in selected NP scenarios. In the last column, the zeroes correspond to the 
scenarios indicated in brackets. All values are given in GeV^. In the calculation of the 
large recoil zeroes, we use the Wilson coefficients given in Table [H and rub = 4.68 GeV. 
The zeroes of ApB coincide exactly with those of P2, except in the presence of scalars. 



As mentioned before, these zeroes cannot produce any singular behavior in the co- 
efficients Jj. The only potential singular points are the zero of P2 and the zero of the 
combination 3 + Pe + 3M2. However, 3 + + 3M2 is always strictly positive, and the zero 
of P2 coincides exactly with the zero of the forward-backward asymmetry, and cancels 
out. This is true in the absence of scalars; in general, the relevant parameter is Xj which 
is always well defined, as discussed in Section [51 

The observables M2, 5*1 and S2 are affected by scalar Wilson coefficients (see Fig. [2]) 
- we recall that Cs effects in P5 and Pq are negligible-. Moreover, M2 plays a special role 
since it is the only observable in i? — t- K*ii sensitive to Cp. As can be seen in Fig. O the 
most promising place to look for scalar effects is in the observable 5*1, when integrated in 
the full region. This is due to the fact that for fixed Cg\ the curves are always positive 
(if Cs > Cg) or negative (if Cs < Cg), while in the absence of scalars Si is strictly zero 
in the full range. Integrating over has a clear experimental advantage from the point 
of view of statistics. However, we would need to discriminate values for the integrated 
observable below the level of | / dq'^Si\ ~ 0.12 GeV^ in order to improve the Bg fJ^'^f-i 
bound. 

In the case of M2 , considering the fully integrated region also leads to better sensi- 
tivity to C^p , since M2 > for all g^. However, the sensitivity in this case should be better 
than A(J dq^M2) ~ 0.03 GeV^, before the current bounds on the difference Cp — Cp can 
be probed. 

Finally, for the sake of completeness, we show in the third row of Fig. |2]the observables 
dV /dq^ and A-pB in the SM and in Scenarios A, B, C and D. It becomes manifest that these 
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Observable 


Wilson Coefficients 


Pi 




P2 


C7, C'-j 


P3 




Pi 




P5 


C'lQ, [Cg] 


Pe 


Im(cW,cW,cS) 


Ml 




M2 


[Cp — C'p] 


Si 


Cs - C's 


S2 


[Cs - C's] 



Table 6: Main contributions to the observables from NP Wilson Coefficients. The listed 
WC's produce strong deviations from the SM (consistent with all other bounds). For 
those hsted in brackets, the effect is milder. 

observables are affected by larger hadronic uncertainties, and show a milder sensitivity 
to NP effects. Although the forward-backward asymmetry seems to suffer a significant 
enhancement in Scenario A. 2, this scenario can be much more effectively probed by the 
observables Pi or P4. 

A summary of the NP sensitivity of each observable can be found in Table |6l For each 
observable, we list the Wilson coefficients whose NP contributions affect most strongly the 
values of the given observable. We also present within brackets those Wilson coefficients 
whose effect is moderate. 

8 Summary of results 

The angular distribution of the four body decay Bfi — j- )■ K7r)i~^i~ can be studied 

experimentally by doing a fit to the coefficients Ji{q^) of the distribution, defined custom- 
arily as in Eq. ([1]). The contact with theory is given by the expressions of the coefficients 
Jj in terms of the transversity amplitudes, as shown in Eq. ([3]) [where fii is defined in 
Eq. (jl])]. In general, these amphtudes are A^\^q-, At and As- 

However, depending on the case (if the masses of the leptons are negligible, as for 
example if = e^, or if there are no NP contributions from scalar operators), these 
coefficients are not independent. In such cases, an independent fit to all the coefficients 
can be problematic. Moreover, since such correlations contain physical information, it is 
interesting not only to have them identified, but to take profit of them. 

On the other hand, the coefficients J, are not the best observables to consider from 
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the theory point of view because they suffer from large hadronic uncertainties. This has 
been noticed before and many theoretically clean observables have been devised in the 
literature. However, not all the observables that one can construct from the transversity 
amplitudes can be extracted from the angular distribution if they violate some "symmetry 
properties" . 

With the development of a formalism based on these "symmetries" , together with the 
above considerations, in this paper we have constructed a complete and efficient set of 
observables engineered to extract the maximum information from the angular distribution: 

1. In the most general case, the chosen basis of observables is composed by the FFD 
observables dV/dq^ and A^-q [Eqs. (fTT]).( l46|) ]. and the FFI observables -Pi,2, 3,4,5,6 
[Eqs. (HI- (HID], Mi,2 [Eqs. (I27D,(I28D] and ^1,2 [Eqs. (|42D,(l43D]. The angular dis- 
tribution in terms of the observables is given by Eqs. ( j3Tl) - (l40l) . with x defined in 
Eq. (HTll . The uniangular distributions can be found in Eqs. (!52l)- (l5^ . 

2. The reduction to the As = case is obtained by setting Si = and x = ^fb/-P2- 
The vanishing of Si leads to two dependencies among the J's: Jqc = and the 
relationship of Eq. ([5]). In fact, 5*1 and S2 measure the breaking of these relations 
by scalar effects. 

3. The reduction to the massless case {nii = 0) is obtained by setting — 1 and 
Mi = 0. This leads to two further relationships between the J's: Jis = 3J2s and 
Jic = —J 2c- In fact. Ml and M2 measure the breaking of these relations by mass 
effects. 

The NP sensitivity analysis shows that these observables are quite sensitive to com- 
plementary NP effects. This can be observed in Figs. [Tp] and in Table [61 It is almost 
a certainty that future analyses of LHC data by the LHCb collaboration will be putting 
serious constraints on NP physics by studying these observables, or else discrepancies 
with respect to our SM predictions will be made manifest and constitute part of the first 
studies of true physics beyond the Standard Model. 
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A Symmetry Formalism 



In this appendix we complete the symmetry approach to the angular distribution that was 
originally presented for the massless case in Ref . [7] . We present two different formalisms 
to describe the distribution. The first formalism, constructed using unitary matrices 
and two-component complex vectors, will be appropriate to describe both the massless 
and massive cases. However, in order to introduce the scalar contributions a more general 
formalism is required. This second more powerful formalism introduces, instead, ortogonal 
matrices and four- component vectors and is valid for all cases. 

We follow here a bottom-up approach, from the simplest (massless) case to the gen- 
eral case (massive with scalars). We also recall the solution of the system, expressing 
transversity amplitudes in terms of J's, in the massless case (see Ref. [7]), while the full 
solution to the system in the general case will be presented elsewhere [35] . 

The importance of determining these symmetries is mainly twofold. On the one hand, 
from the experimental point of view, the symmetries allow to identify all correlations be- 
tween the coefficients of the distribution that may affect the stability of the fit; but they 
are also helpful to determine which amplitudes can be consistently put to zero, in order to 
simplify the system and consequently the fit. On the other hand, they provide you with an 
alternative procedure to construct observables directly in terms of the transversity ampli- 
tudes: verifying that they are invariant in the first place, and afterwards, translating their 
expression in terms of transversity amplitudes to an expression in terms of the measured 
coefficients Ji of the distribution (an example of this procedure was the observable 
designed in [7]). 



A.l Symmetries of the massless distribution 

In this section we review the symmetry formalism for the massless angular distribution, 
as presented originally in Ref. [7]. 

The six complex amplitudes present in this case can be arranged into three complex 
vectors: ^ 

^11 = f fil ' ^^={ . ^0 = . (58) 



All the coefficients Jj can be expressed in terms of the products n\ nf 

Jis = ^ (|n_L|^ + |n|||^) , Jic=\no\^, -/2s = ^ + jny , 

^2c = -|no|^ = ^ (I'^iP - I'^llH , -''4 = -^Re(raJnii) , 

J5 = -\/2 Re(?7,o n_L) , Jqs = 2 Re(?T,|j_ ny) , J7 = —V2\m.{n\n\\) , 

Js = ---^Im(n|)n_L) , Jg = -Im(n|L riy) , Jec = . (59) 
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A symmetry of the angular distribution will therefore be a unitary transformation U 
acting in the same way on rio, n\\ and n±, that is: rii — > Uui. Such a symmetry has four 
independent parameters, and can be written as: 



Urii 



o1-<PL 











-i<t>R 



COS 6 — sin 6 
sin 6 cos 6 



cosh id — sinh i6 
- sinh i6 cosh 



Hi 



(60) 



Of course, other parametrizations are possible, but we prefer to keep this one to make an 
easy contact with the generalization to the massive case and the notation introduced in 
Ref . [7] . The matrix U defines the four symmetries of the massless angular distribution: 
two global phase transformations {(pL and 4>r), a rotation 6 among the real and imaginary 
components of the amplitudes independently and another rotation 6 that mixes real and 
imaginary components of the transversity amplitudes. 



A. 1.1 Solution to the massless distribution 

We can now use these symmetries to reduce the number of theoretical parameters and 
solve for the transversity amplitudes in terms of the coefficients J's. It is instructive to 
use only three out of the four symmetries and see how the extra freedom related to the 
fourth symmetry arises from the equations. This extra freedom gives rise to the non-linear 
relation between the J's given in Eq. (I5j). 

Using the symmetries we choose to fix the following amplitudes to zero: = and 
ImAji^ = 0. We achieve this configuration easily using the rotation phases (pL and (pR to 
set the phases of Af^ and Af^ to zero. Then a rotation by an angle 6 given by 



tan^ 



ReAf 
RelF 



(61) 



will also set the modulus of A^^ to zero. 

With these simplifications, rewriting the products nl rij in this configuration, and 
taking into account Eqs. fl59|) . one gets immediately: 



A^ = 



A 



R 



A 



R 



Jpis — 2z Jo 



A 



R 



\/lJ2s — Js 

2 J4 — iJ-j 



I^Lhs - J3 V4J2S - 2 J3 

The remaining equations involving the last two amplitudes {^A\_ and ) lead to [7] 

_ 2(2 J2, - J3)(J5 + lih) - (2J4 + ?J7)(^6s - 2?J9) 



(62) 



4J| 



■Hs ~ 4J|a/2 Jic(2J2s — J3 



4 J| - J| 



(63) 



where 0q and 



\ are the phases associated to the amplitudes A^ and A\. The relation 
between the J, coefficients (Eq. ([5])) arises naturally from imposing in Eq. (!63|) that the 
modulus of this phase difference should be one. Notice also that the freedom to choose 
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one of the two phases in Eq. (1631) to be zero is somehow related to the freedom associated 
to the last unused symmetry transformation 6. The choice 0;^ = fixes the last two 
amplitudes to 

2(2 J2. - + 2zJ8) - (2J4 + iJi^hs - 2zJ9) 

V4J2S - 2J3 V16J|, - 4J| - Jl - 4J| 

The solution in any other configuration can be obtained by applying the symmetry trans- 
formation in Eq. ( I60l) . Any observable constructed from the transversity amplitudes, can 
be expressed in terms of the coefficients Jj using Eqs. (!62|) and The condition that the 
observable is invariant under the symmetries of the angular distribution, guarantees that 
the configuration used to derive these equations gives the same result for the observables 
as any other configuration, and the result is unique. 



A. 2 Symmetries of the massive distribution 

In the massive case the balance equation between theoretical and experimental degrees of 
freedom [7] 

nexp = nj -nd = 2nA - ns (65) 

is fulfilled with rij = 12, Ud = 2, ua = 7 and ns = 4, where Ud is the number of 
relationships among the Jj, as explained in Section |2l When the masses are switched 
on to account for lepton mass corrections in Bd — ?■ /i*°(— )■ K7r)i^i~, all the massless 
symmetries described in the previous section are broken by the mass terms and have to 
be redefined. The four symmetries in the massive case are then: 



A common global phase transformation for both left and right components 



n'i = Uo{(f))ni 



e'* 
e"**^ 



(66) 



for i = II, -L, 0. 

An independent new global phase transformation for the At amplitude: A[ = e^'^^At. 

Two rotations Ui{6) and 1/2(0) between the real and imaginary components of the 
transversity amplitudes Af'^: 

n[ = Ui{9)ni, n[ = U2{0)ni, 

{i =_L, 11,0) with a similar structure to those in the massless case, but including 
some important differences to be discussed below. 
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In order to find the explicit form of the last two symmetries, Ui{d) and U2{0), it is 
helpful to analyze their infinitesimal form. This can be obtained following the approach 
described in Ref. [7]. Let us focus on Ui{6) and require that the infinitessimal transfor- 
mation of the amplitudes is a symmetry of the distribution. This leads to the following 
form of the infinitessimal transformation associated to a rotation with angle ^ ~ e, 



A 



Li 



Li 



A 



^0 



A 
A 



L ^ aR* 



A 



Ri 

Af 
A^ 
\A 



aL _ ^ aR* 



Af - eAf* - leAfk , 
+ eAf * - ieAfk , 
Aq -\- (^Aq ^^^0 ^ ' 

.l2 





/ 12 
\L2\ 



At\'' - 2e [Re(A^O - Re{Af) + klm{Ai^' A^)] , 



(67) 



where k = [(Re(Af 2) - Re{A^^)) - (|| o±)]/[Im(Aj'Af ) + (|| Eq. (EZD together 

with some important observations will guide us in the construction of the corresponding 
continuous transformation, namely: 

I. The structure of the distribution and the absence of lepton masses in the infinites- 
simal transformation informs us that the symmetries of the massive distribution 
should be also symmetries of the massless case; consequently, the form of the rota- 
tion matrices should respect the form of the transformations in the massless case. 

II. The infinitessimal form shows that while the left components transform linearly, the 
right components transform non-linearly. Moreover, in the limit of /c — )■ 0, the linear 
transformations of the massless case are recovered. 

III. All Ji, except for Jis and Jic, are invariant under this infinitessimal transformation 
independently of the explicit form of k. This last remark is, indeed, connected to 
point I. 

These considerations taken together imply that the continuous symmetry transformations 
Ui{6) and U2{0) should take the form: 



Ui{9) 



cos 6 — sin 6 



cosh iO — sinh iO 
.g-i<5(9) sinh 26* e-^^(^) coshie^ 



(6^ 



The last step to determine these rotation matrices completely is to obtain the phases 6{6) 
and S{6). These are non-linear functions of the transversity amplitudes and the angles 
6 and 6, respectively. In a certain sense, the goal of these non-linear pieces is to cure 
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the breaking of the massless symmetry by the massive terms, while respecting the basic 
structure of the massless symmetry. Imposing that the mass term in Jig 

Re {AiAf* + A^Af) (69) 

should be invariant under the symmetry transformation, we obtain that sin 5{6) and 
cos S{6) are just the result of a rotation of a unitary vector (sinw, cosv) whose first com- 
ponent (sin v) is indeed proportional to the mass term. This rotation is 



sin5(6')\ / cosM — sinw 
cos 6 {9) J \ sinu cosu 




(70) 



where sinf = xi/ \/h{6), cosv = riy^^Jl — x\/h{6) and r]y^ is the sign of the function 
Hi. The sign function has been introduced to ensure that 5(0) = 0, which implies that 
the transformation matrix becomes the identity matrix^] for 6* —t- 0. The rotation (n) is 
defined by 

?/i cos 26' + 2/2 sin 2^ Xi cos 26^ + 0:2 sin 20 
cosM = , , smM= -== , (71) 

where 

xi = Re(4Af ) + (II , X2 = l [MAi') - Re(Af 2)] - (|| , 

y, = Im(4Af ) + (II ^±) , y, = ^ [imiAf') + lra{Af')] - (|| ^±) , 

h{e) = {xi cos 29 + X2 sin 2ef + [yi cos 29 + y2 sin 2df . (72) 

Xi is precisely the mass term of Ju. Being an invariant, it can be expressed in terms of 
the coefficients Ji of the distribution: 



_ (2 + /3J) / Ji. .hs . 

""'-^.^^[JVM ^^^^ 



Also the non-linear parameter k in Eq.f l67|l can be rewritten as /c = 2x2/ yi- Notice that 
due to the non-linear form of the transformation 

f/l(^l)-f/l(^2)7^f/l(^l+^2). 

This does not pose a problem since one can easily concatenate transformations, one after 
the other, always evaluating the corresponding 5{9i) for each transformation. A final 
important remark is that the requirement of a/1 — x\/h{6) to be real imposes a restriction 



■''Notice, however, that the cosu defined as cosu — —rjy^ ^Jl — x\/ h{8) is also a solution, even if not 
connected to the identity in the limit — > 0. 
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for the range of validity of the transformation around ^ = 0, i.e., there is a maximum and 
a minimum allowed value for 6*, given by the condition 



Xi 



< 1 . 



(74) 



Inside this range of validity 6{6) is simply given by 6{6) = v — u. 

To complete the transformation of all amplitudes under Ui{6) we still need to find the 
transformation of the amplitude At. This is obtained by imposing the invariance of the 
mass term in Ji^ 



\At\' + 2(Re{A^Al^* 



Re{Af;'Ar: 



(75) 



where the primed Aq'^ can be obtained easily from Uq. 

Following exactly the same steps for the symmetry transformation U2{0) we can also 
identify the corresponding 6 (9) by the rotation 



sin5(^) 
cos 6(6) 



cos n — sin n \ / sin v 
sintt cosM / \cosv 



(76) 



with 



where sin w = xi/ y h{0), cosv = rjy^d 1 — x\/h{6), rjy^ is the sign of the function yi and 



cos u 



iji cos 26 + 2/2 sin 26* 



smu 



Xi cos 29 + X2 sin 26 



(77) 



Xi = Xi 



X2 



[Im(y4p) - Im(v4 



R2 



m =yi, y2 = -[-Re(4')-Re(A 



R2 



h(e) = (xi COS 29 + X2 sin 29f + {yi cos 29 + y2 sin 29^ . 



(7J 



The same discussion about the range of validity can be applied to U2{9) just by substi- 
tuting in Eq. (I71D h{9) — h{9) and now 5{9) = v — u. \A[\'^ can be computed 
from Eq. (!75|) but using the corresponding primed Aq'^ amplitudes under the U2{9) 
transformation. 



A. 3 Symmetries in the presence of scalar contributions 

The last case to discuss here is the symmetry structure of the angular distribution in 
presence of scalar contributions to the decay channel — )■ K*^{-^ Kti)£^1~ . 

The infinitessimal transformation, too long to write it explicitly here, provides the 
following important information: 
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• A± and A\\ transform exactly as in the massive case. 

• Aq and As get mixed in the transformation and, contrary to the massive case, an 
exphcit dependence on the lepton mass appears in the symmetry transformation. 
This has an important consequence for the construction of observables: while the 
lepton mass terms in Jis is invariant by itself, the mass terms in J5 and J7 are not. 

Another fundamental difference with the previous cases, is that the use of the compact 
two-component complex vector rii is no longer possible for Aq and As- Consequently, we 
will introduce a new formalism in terms of four- component vectors 



ImA^ 
ReAf 
\lmA^ J 



v± 



( ReA^ \ 
ImA^ 
-ReA^ 
y -\vclA\ ) 



( ReA^ \ 
ImA^ 
ReA^ 

l^ImA^ j 



vs 



( ReAs \ 
ImAs 





(79) 



and two 4x4 matrices 



C 



I2 I2 
I2 I2 



7 



-io2 

i(T2 



and will describe the angular distribution in terms of them. Notice that 7^ 



7-7 



I4. We remind here the explicit form of the Pauli a matrices: 



1 

1 



0-2 



-2 

1 



1 
-1 



(80) 
—7 and 

(81) 



The matrix C is needed to symmetrize the vector v*s^ which appears in the angular 
distribution as vsc = C ■ vs- In terms of the four vectors in Eq. (179|l and the two 
matrices in Eq. (l80|) one can rewrite the coefficients Jj including scalars as 



Jls 
Jlc 

Js 



(2 + /3|) 



vl-vl + v\i ■vli) + 2^{v\i ■ {C - I4) ■ t^i - t'l ■ {C - I4) ■ vl] 



, m 



/3| 



vo ■ iTq + 4^{vo ■ {C - I4) ■ ^To + I Ar) + ^{vsc ■ vsc) 

^(wl ■ Wl + V\\ ■V\(), J2c = -Pl{Vo ■ Vo) , 



2 



^ (wl ■ v± — m ■ v\\) 
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me 



vl-'y ■ Vsc) , 



'V2 



/3]{vl- -f- Vo) , Jg = ■ 7 ■ wl) 



The angular distribution exhibits four symmetries: 
• A global phase transformatioiil^ for 

v'j = V(o)(0)w;-, 

with j = II, J-, 0, SC, and where 



(0) 



cos </) — sin </) 

sin0 cos</) 

cos — sin 

sine/) cos 6 



(82) 



(83) 



U) 



• An independent phase transformation for the At amplitude: A[ = e^'^^At . 

• And the same two rotations of the massive case. However, there is an important 
difference between the transformation properties of A^'"^, A^^'^ on one side and the 
transformation properties of v4q'^, ^5 on the other, that will be detailed in the 
following. 

On the one hand, A^^^ and A^^'^ transform exactly as in the massive case 



V{i){6)vj and 



Vi2){e)^ 



i5) 



for j =±, II and the matrices V(i) (6) and V(2) (6') are simply the mapping of the 2x2 matrices 
Ui{6) and U2{0) in the 4-d formalism and are defined by 



Vii){9) 



cos 6 — sin 6 

cos 9 sin 6 

cos 6 sin 6 sin 6 sin 6 cos 6 cos 6 — sin 6 cos 6 

sin S sin 6 — cos 6 sin 6 sin 6 cos 6 cos 6 cos 6 



(86) 



^Incidentally notice that the matrix 7 can be interpreted also as a phase transformation of 7r/2 for 
the L components and -7r/2 for the R components. 
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and 



COS 9 — sin 9 

cos 9 — sin 9 
■ sin 5 sin 9 cos 5 sin cos 5 cos — sin 5 cos ^ 
cos (5 sin ^ sin (5 sin ^ sin 5 cos ^ cos 5 cos 6' 



(87) 



The non-linear structures 5{9) and 5{9) in V(i){9) and V(2)(6') are the same as in Eq. ( ITOl) 
and Eq. f l76|) respectively. In this new formalism, the expression of Xj, Ui, Xi and iji can 
be rewritten in a more compact way: 



(89) 
(90) 
(91) 



and 



with 



Xi 


= {n ■ 




-I4) 


■ V\\ - 


- vl 


{C-I,)-v^) /2 




= in ■ 


Cx2 


■ V\\ - 


- v\_ 


■ Cx2 


■v±) /2, 


yi 


= {n- 


Cyl 


■ V\\ - 


- ^'1 


■ Cyl 


/2, 


2/2 


= (-ii- 


Cy2 


■ V\\ - 


- ^'1 


■ Cy2 


■ VA.) /2 




Xi = 


Xi , 




yi = 


yi , 






X2 = 


(-11 


Cx2 


■ V\\ - 


- vl 


Cx2 ■ VA_) /2 , 




m = 


(^'11 


Cy2 


■ V\\ - 


- ifl ■ 


Cy2 ■ V±) /2 , 



(92) 
(93) 
(94) 



C 



x2 




c 




C 



3/2 




c 



x2 




Cfi 



y2 




(95) 



Li Ft 

On the other hand, in the scalar case the transformation of the amplitudes ' is different 
from the massive case, but shares the same structure than the transformation of As: 



where 



-03 
\vofj 



^'sc = (vs^^^ + Vsc 



(96) 



VS2 

\vs?J 



(97) 
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with j = 1,2. These vectors contain the non-hnear part of the transformation associated 
to Aq'^ and As and are functions of 6 (for j = 1), 9 (for j = 2), all amplitudes and the 
lepton mass me. 

The last remaining point is to determine Vq^-''* and vs^-'\ They are fixed by requiring 
the invariance of the angular distribution under the Vi(^) and ¥2(6) transformations. The 
set of eight equations required to obtain the components of Vq^-''^ and vs^^"* arise from the 
coefficients J4, J5, J7 and Js after imposing the invariance of the angular distribution 
under the transformations: 

v\\ ■ Vq^^^ = , w"! ■ 7 ■ Vq^^^ = , 
vl ■ vq^^'^ - ^v\\ ■ vs^^^ = , ■ 7 ■ vo^^^ - -^vl ■ 7 ■ vs^^^ = , (98) 



together with the ones derived from and Jgc 

y^U) . y^U) + 2vo^^^ ■Vo = , vo'-'^ ■ vs^^^ + vo^^^ ■ vsc + Vo ■ vs^^^ = , (99) 

and two more equations to impose that the first component of the vector vsc must equal 
the third while the second component should equal the fourth^ 

Jti) . (.r^O) + y-^) = , . (tT^O) + Vsc) = , (100) 

where 

(^^^ = - ^(i)3,l, ^0)1,2 - V(j)3,2, V(j)i,3 - V(j)3,3, Vq-) 1,4 - V(j)3,4) , 

e^^^ = {V{j)2,l - V(j)4,l, V(j)2,2 - V(j)4,2, V(j)2,3 " V(j)4,3, V(j)2,4 - Vl^j) 4,a) ■ (101) 

The system can be solved as follows: 

I. Obtain ■Uq*'"'^ in terms of Vs^^^ from the first set of Eqs. fl98l) . The components of the 
vector Vq^^^ are then: 

vq[^^ = mS^'^ReX^'^ + m^^'^ImF^'^ , wq?^ = m[^hmX^'^ - m^^'^ReF^'^ , 

v^^p = -m['^ReX^'^ - mi'hmY^^^ , Voa^ = -m?lmX^^^ + m^^^ReY^'^ , (102) 

where 

X^''' = AjjAfAj + \A\\^Al, 

yaf, ^ ^o^^*b^b^ ^ lA'^I^Ajj , (103) 

with 



^Therefore vsc' = {ReA's,lmA's,RcA'g,lmA's). 
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and w = \A^\^\Af\^ + \Af\^\A^\^ + 2Re{ A^* A^A^Af). 



II. Use Eq.f llOOp to express the components ^53"''' and vs^' in terms of vs^^' and 1^52 



(i) 



(i) 



(i). 



0? 





- dj'^ei^') r 


"4 '^s 


_ jU)pU) _ 




- 4'^4'^ r 


"4 '^s 


_ ^(i)„(i) . 

"3 ^^4 



+ 



d. 



1 Co 



"4 ^^3 



"3 ^^4 



i;5^^') + ImA5 



"3 ^^4 



ImA.c 



-ReA^ 

- ImAs 
(105) 



III. Determine vsi^ and ■ Using Eqs. fl99l) together with Eq. fll02p we find the last 
two equations. The first of Eqs. f l99|) gives 



(i)2 
1 



yL,ij|2^0)2 ^ 2mi^^m^^^Im(X^'^*F^'«)] + (L ^ i?) 



-2Re(A^*X^'^)mi(j) - 2Im(A^*F^'^^ 



(i)i 



(L ^ i?) 



(106) 



where are functions of vs^^"* as given in Eq. fll04p . 
The second of Eqs. (!99|) gives rise to 

[mS^)Re(^JX^'^) + mi'hmiA*sYL,R)] - [L ^ R) 



(107) 



-^^[^'^(ReA^ + mi^^ReX^'^ + m^^^ImF^-^^ 



-vs2\^^^h + m^^ImX^'^ - m^^'^ReF^'^) 
-vsfi^^A^ - mS^'^ReX^'^ - m^^'^ImF^'^) 
-vs'i\^^A^ - mfhmX^^^ + m^^^ReF-^'-^) 



IV. Substituting Eqs. (fT04ll and Eqs. (|T05ll in Eq. (fTOHjl and Eq. (ITOTll we end up with 
a system of two coupled quadratic equations which are function of vsi'' and ^52"'^ 
This system can be solved numerically and typically provides two complex solutions 



(to be discarded) and two real ones. The real solutions for vsi^ and vsi^ once 



(i) 



inserted in Eq. ( 196|) . generate two sets of transformed amplitudes v'q and v'g^ that 
leave the angular distribution invariant. One of them is connected to the identity 
whereas the other is not, exactly as it occurred with 5{9) and 5{9). This completes 
the definition of the symmetry transformation of the Aq ' and As amplitudes. 



34 



B Large recoil limit expressions 

In this appendix we present the expressions of the observables Pi and Mi in the large 
recoil limit. These are useful to study qualitative properties of the observables as well as 
for rough quantitative estimates. 



|n_L 
Re(?7,^n||) 
Im(?T,^n| 
Re(?T,Q?T,| 



Im(nQn||) 
Re(?T,Qnj_) 



|?7,_|_P — |?7,|| 



I |2 I I |2 

\n\\\ + n_L 



Ml 



Mo 



pt 1 2 
'-'lOl 



2d-2 (\p nl |2 



+ Cf + Cf) - {Fs Cf + Cf )!') 
|(FCf + Cf) - (FCf + Cf )!') , 



2f2 /ir I r" |2 



'10 I 



I (FCf + Cf) + (FCf + Cf )l') 



4Q?elRe[(FCf 



Cf )C*o - (FCf + Cf )C 



lOj ' 



-4Q2^1lm[CioC;* + (FCf + Cf )(FCf + Cf )*] 



2Qfg2eiie± (|Cio 

+F(1 + s)Re[(Cf - Cf'){Cf - Cf')*]) 
-2Q?Q2e||aim[F(l - s){Cf - Cf )(Cio - C 



cff _ r>cS/ 1 2 
7 '-'7 I 



2QiQ2^ll^±Re[i{F{l + s)Cf + C 



-.cff\ 



-((F(l - s)C 



off 



7 "T "-9 

Cf) + (F(l + S)Cf + Cf))Cioj, 



7Cff/ 

'7 

7Cff/\\/l^* 1 



(F(l-s)Cf + Cf))Cro 



(108) 
(109) 
(110) 
(111) 
(112) 

(113) 
(114) 

(115) 



Im(4nx) = 2QlQ2^\\^± (2Im[CioCro + F^s CfCf* + CfCf'*] - Im[F((l - s)Cf 



+ (1 + s)Cf')Cf* - F((l + s)Cf + (1 - s)Cf')Cf'*]^ , 

8Qia (RepioC;*] + Re[(FCf + Cf){FCf' + Cf )1) , 
4Q?el (ICior + |C(or + \FCf + cf |2 + |FCf + Cff) 



2mj |Cio 




""-{FCf + Cf 


""-iFCf' + Cf' 


2 


|Cio 


^ + \C[o 


2 + |FCf + Cf 


^ + \FCf' + Cf 


2 



4mi 



where the following short-hand notation has been used 

^ 2mb „ /:r „ 1 



Qi = y2iVmB(l-s), Q2 = 



2\f2mK*^ 



(1 



(116) 

(117) 
(118) 

(119) 
(120) 

(121) 
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